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ABSTRACT 

We have compiled a sample of 234 ultra-steep-spectrum(USS)-selected radio 
sources in order to find high-redshift radio galaxies (HzRGs). The sample is in the 
southern sky at —40° < 6 < —30° which is the overlap region of the 408-MHz Revised 
Molonglo Reference Catalogue, 843-MHz Sydney University Molonglo Sky Survey (the 
MRCR-SUMSS sample) and the 1400-MHz NRAO VLA Sky Survey. This is the third 
in a series of papers on the MRCR-SUMSS sample. Here we present optical spec- 
tra from the ANU 2.3-m telescope, ESO New Tec hnology Telescope and ESO Very 
Large Telescope for 52 of the identifications from iPaper III , yielding redshifts for 36 
galaxies, 13 of which have z > 2. We analyse the K-z distribution and compare 4- 
arcsec-aperture magnitudes with 64-kpc aperture magnitudes in several surveys from 
the literature; the MRCR-SUMSS sample is found to be consistent with models for 
10 11 -10 12 M Q galaxies. Dispersions about the fits in the K-z plot support passive 
evolution of radio galaxy hosts since z > 3. By comparing USS-selected samples in 
the literature, we find that the resultant median redshift of the samples shown is 
not dependent on the flux density distribution or selection frequency of each sample. 
In addition, our finding that the majority of the radio spectral energy distributions 
remain straight over a wide frequency range suggests that a k-correction is not respon- 
sible for the success of USS-selection in identifying high redshift radio galaxies and 
therefore the steep radio spectra may be intrinsic to the source or a product of the 
environment. Two galaxies have been found to have both compact radio structures 
and strong self-absorption in the Lya line, suggesting they are surrounded by a dense 
medium. For the bulk of the sources, spectral line ratios show that photoionisation is 
the primary excitation process. 

Key words: galaxies: active - surveys - galaxies: high-redshift - infrared: galaxies - 
radio continuum: galaxies. 



1 INTRODUCTION 

Hierarchical models of galaxy formation predict that mas- 
sive galaxies formed through many mergers of smaller sys- 
tems, with continuous bursts of star formation. However, 
there is solid evidence that the most massive galaxies in the 

* Based on observations obtained at the ESO Very Large Tele- 
scope (VLT) and New Technology Telescope (NTT) as part of 
program numbers 077.A-0471 and 079.A-0504 and the ANU 2.3 m 
telescope. 

f E-mail: jbryant@physics.usyd.edu.au 



local Universe for med their stars very early in the history 
of the Universe ([Stanford. Eisenhardt fc Dickinson 1 1 19981 : 
iNelson et alll200lt iThomas et alj|2005l h In order to fit cur- 
rent galaxy-formation models, feedback processes have been 
introduced in which the energy input from a central mas- 
sive black hole quenches the star formation by heating the 
interstellar gas l|Springel et al.l 120051 : |Pe Lucia et al.l 120051 : 
iBower et al.ll2006h . Such modifications appear to overcome 
many of the shortcomings of previous hierarchical models, 
but although the quenching of star formation is attributed 
to the onset of a galactic wind driven by ionising radiation 
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and/or mechanical energy from the radio jets, the physi- 
cal mechanisms for this are poorly understood. In contrast, 
there is observational evidence to suggest that the radio jets 
associated with some supermassive black holes might trig- 
ger, rat her than suppress, star formation in high-redshift 
galaxies (lReeslll989l ; lBicknell et al 1l2000l ; iKlamer et al]|2004 
ICroft et alj2006h . There are also plausible arguments to sug- 
gest that both processes operate, with radio jets first trig- 
gering star formation and then later driving a galactic wind 
The interplay between mergers and feedback as- 
sociated with active galactic nuclei (AGN) is likely to affect 
the evolution of all massive galaxies, not just the extrem e 
bright end of the luminosity function (|Croton et alj [2006). 
However, it is the most massive galaxies at high redshift that 
provide the most stringent constraints on these processes. 



The star formation rate density in the early Uni- 
verse is dominated by the mo st massive galaxie s in 
the sample volume (see fig. 2 of iFeulner et al] 120051 ). In 
iRocca-Volmerange et al] (|2004l . fig. 5), radio galaxies have 
been identified to be more massive than optically-selected 
galaxies at each redshift, with stellar masses of ~ 10 12 M Q , 
which is 5 — 10 times more massive than the largest galax- 
ies in most optical surveys. The star formation rate den- 
sity of the Universe has been shown to p e ak at 2 < z < 
3.5 jLillv et al] Il996l; iMadau et al] Il996l ; iHopkinsl 12004 
iHopkins fc Beacoml 20061 . fig. 2) and since radio galaxies are 
the most massive and dominate the star formation rate at 
these epochs, they are crucial to understanding early galaxy 
evolution. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 Target selection and radio and _K-band 
imaging 

The target selection, radio observations and K-bwad imag- 
ing were discussed extensively in Papers I and II. A very 
brief summary is given here. 

The MRCR-SUMSS sample was formed by cross- 
matching the revised 408-MHz Molonglo Reference Cata- 
logue (MRCR), the 843-M Hz SUMSS catalogu e and the 
1400-MHz NVSS catalogue (|Condon et al]|l998h in the re- 
gion where all three catalogues overlap (declination —30° to 
—40°). We then selected sources with a spectral indejQ be- 
tween 408 and 843 MHz, atol < -1.0 and Galactic latitude 
\b\ > 10°. We further eliminated any sources that had an- 
other NVSS source within 100 arcsec to minimize the effects 
of source confusion. 

Follow-up radio images were obtained with the Aus- 
tralia Telescope Compact Array (ATCA) at 1384 and 
2368 MHz for the full sample, to pinpoint the isT-band iden- 
tificatio n and in vestigate the radio structure. These are de- 
tailed in lPaper J . Twenty-nine of the MRCR-SUMSS sources 
were also o bserved a t 4800 and 8640 MHz; those data are 
presented in iPaper III 

7f-band (2.2/^m) observations were made with the 
3.9-m Anglo-Australian T elescope's (AAT) I RIS2 detector 
ilGillingham fc JonesfeoOOT ) and with PANIC l|Martini et all 
l2004h on the 6.5-m Magellan Baade telescope at Las Cam- 
panas Observatory. F ull details of the A"-band imaging and 
reduction are given in iPaperlJ . 



We have conducted a survey of ultra-steep-spectrum 
(USS) radio sources in the Southern Hemisphere, designed 
to find high-redshift (z > 2) radio galaxies so as to carry out 
a detailed quantitative study of their environments. This is 
the first USS survey in the South to use large, deep low- 
frequency catalogues to mine a large area of sky, and it be- 
came possible with the completion of the S ydney University 
Molonglo Sky Surv ey at 843 MHz (SUMSS; lBock et al.lll999l; 
iMauch et all 120031 ') and the reanalysis of the 408-MHz Mo- 
longlo Cross survey ( Large et al.lll98lT ) to give the more sen- 
sitive revised Molonglo Reference Catalogue (MRCR; Craw- 
ford, private communication). 



The details of the sample selection, and the 1384- and 
2368-MHz r adio imaging of the MRC R-SUMSS sample were 
presented in lBroderick et ail (120071 . hereafter Paper I). The 
A'-band imaging and identification of the co unterparts to 
175 o f the radio sources were presented in iBrvant et al] 
(2009, hereafter Paper II). The current paper includes the 
followup spectroscopy for the sources in IPaper III and a com- 
plete discussion of the sample using the radio, Jf-band, red- 
shift and spectral line information. A following paper will 
discuss the environments of the z > 2 galaxies and what 
that can tell us about the early evolution of massive galax- 
ies. 



2.2 Spectroscopy 

Spectra for our sources were obtained using three spectro- 
graph/telescope combinations. Sources with o ptical identifi- 
catio ns from the SuperCOSMOS Sky Survey l|Hamblv et al] 
2001) were observed using the 2.3-m telescope of the Aus- 
tralian National University (ANU) at Siding Spring Obser- 
vatory, NSW. Identifications with K <19 were observed with 
the 3.6m ESO New Technology Telescope (NTT), and the 
faintest targets (K > 19) were observed using the 8 m ESO 
Very Large Telescope (VLT). For both the NTT and VLT 
observations, to acquire the faint targets into the slit, we 
used blind offsets from a nearby (<1 arcmin) star in the K- 
band images. After a first exposure, we performed a quick 
data reduction. If the redshift was obvious, the next ex- 
posure was aborted to save observing time. The individual 
exposures were offset by 10 arcsec along the slit to allow sub- 
traction of the fringing in the detector at red wavelengths. 
The journal of observations is given in Table [T] 



2.2.1 ANU 2.3m 

We used the Dual-Bea m Spectrograph (DBS; 

iRodgers. Conrov fc Bloxham I I1988T ) on the ANU's 2.3- 
m telescope on 2007 April 10-14. A plane mirror was used 
to direct all the light into the blue arm of the spectrograph 
in order to obtain complete wavelength coverage over the 



We have adopted a flat, A cold dark matter cosmology 
with H = 71 km s" 1 Mpc" 1 , f2 M = 0.27 and f2 A = 0.73. 



1 Radio spectral index a is related to S„, the flux density at 
frequency u, by S 
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Table 1. Log of observations. 



Date 


Telescope 


Instrument 


Seeing 








(arcsec) 


2006 June 22-23 


VLT 


FORS2 


0.5-1.3 


2006 July 27 


NTT 


EMMI 


1.0-1.4 


2007 April 10-14 


ANU 2.3 m 


DBS 


1.5-3.0 


2007 July 10-12 


NTT 


EMMI 


0.6-1.0 



range 3500-11000 A. The CCD is an E2V 2148 x 562 pixel 
detector with a spatial scale of 0.91 arcsec pixel -1 . We used 
the 158R grating, providing a dispersion of 4 A pixel -1 and 
a spectral resolution of 8.9 A. Severe fringing at the red 
end restricted the usable wavelength range to 3500-8000 A; 
second order contamination was minimal. Conditions were 
variable, but mostly non-photometric, with poor seeing. 

2.2.2 ESO/NTT 

We used the ESQ Multi-Mode Instrument (EMMI; 
iDekker. Delabre fc D'Odoricol Il986h on the NTT on 2006 
July 25-27 and 2007 July 10-12. The conditions in 2006 were 
poor with the first two nights lost to weather; the third was 
non- photometric with 1.0-1.4-arcsec seeing. In 2007 the con- 
ditions were better, with seeing between 0.6 and 1.0-arcsec 
FWHM. We used grism #2 and a 1.5- or 1-arcsec slit; the 
latter was used in the better seeing conditions. The wave- 
length coverage was 3700-9700 A, with a dispersion of 3.5 
A pixel -1 , a spectral resolution of 10.4 A FWHM and the 
pixel size was 0.33 arcsec. To minimise the effects of differ- 
ential atmospheric refraction, we observed the targets with 
the slit oriented at the parallactic position angle. 

2.2.3 ESO/VLT 

We used the FOcal R educer and Spectrograph (FORS; 
lAppenzeller et all |l998) in visitor mode on the VLT on 
2006 June 22-23. We used FORS2 on UT1 (Antu) with 
the 1501 grism and 1-arcsec slit, providing a dispersion of 
6.7 A pixel -1 and a spectral resolution of 21.4 A. The wave- 
length coverage was 3600-8500 A. As the FORS instruments 
have a linear atmospheric dispersion corrector, we did not 
observe at the parallactic angle, and instead the slit position 
angle was chosen to avoid bright stars. 

2.2.4 Data reduction 

Data reduction was performed with the iraf software suite, 
using standard procedures. We removed the bias and pixel- 
to-pixel gain variations from each frame, and then removed 
the cosmic rays using the iraf task SZAP. We extracted spec- 
tra using the iraf task APALL, extracting only the central 
few rows (~ 1 arcsec) to maximise the signal-to-noise ra- 
tio of the spectrum. This optimises our chance of finding a 
redshift, at the cost of not including all the flux from the ob- 
ject, particularly from any extended line emission. The same 
extraction aperture was used to extract a calibration spec- 
trum from the arc lamp, which was then used to calibrate 
the one-dimensional spectrum in wavelength. 

The spectra were flux-calibrated by comparing with the 



spectrum of a spectrophotometric standard taken on the 
same night. No correction was made for Galactic reddening. 



3 RESULTS 

Table [2] lists all the MRCR-SUMSS sources for which we 
have taken spectra. The A"-band images with radi o con- 
tours for these objects were presented in IPaperlj . Some 
sources met our selecti on criteria, but were al so part of the 
SUMSS-NVSS sample l|De Breuck et al.ll2004 Paper I), and 
were therefore not reobserved. For those so urces, the if-band 
magni tudes and positions given are from |Pe Breuck et al] 
(|2004h . 

The columns in Tableware: 

(1) NVSS source name. 

(2) if-band magnitude in a 4-arcsec-diameter aperture. 
Some sources are obscured by foreground stars (see col- 
umn 6) which contaminate the listed magnitude. Sources 
marked SC were visible on the SuperC OSMOS UKR fields 
and no K- band image was obtained fsee lPaper it) . Our anal- 
ysis exclusively used magnitudes in a 4-ar csec aper ture, and 
the justification for this choice is given in IPaperlj . 

(3) A'-band magnitudes in equivalent 64-kpc apertures, 
which were calculated fro m 8-arcsec-a p ertur e magnitudes, 
following the method of lEales et al.l (|l997l ). Only radio 
galaxies (not QSOs) which hav e redshifts and 8-arcsec- 
aperture magnitudes in IPaper ij have been converted. The 
limitatio ns on the se 64-kpc-aperture magnitudes were dis- 
cussed in lPaperTll and are discussed further in Section [47TT2] 

(4) & (5) RA and DEC ( J2000) of the if-band identification. 
One source which had no K-band identification is marked 
" — ". The sources which had SuperCOSMOS identifications 
have the SuperCOSMOS position listed. 

(6) Spectroscopic redshift. The wavelength of each line was 
determined by fitting a Gaussian on top of a continuum. 
Occasionally spurious lines were introduced by the fringing 
in the CCD; we required that features be visible in both 
exposures to be counted as real. The redshift quoted is the 
average of the redshifts derived from individual lines; the 
error is the scatter in this average, with uncertainties due to 
the line fitting and wavelength calibration added in quadra- 
ture. The sources marked 'continuum' had continuum but no 
lines, while those marked 'undetected' showed no continuum 
or line emission. Some objects are obscured by a foreground 
star and are listed as 'obscured'. Objects with emission lines 
broader than > 3000 km s -1 have been classified as QSOs0 
and are marked with a 'Q'. 

(7) Origin of the spectra. On the NTT, EMMI was the spec- 
trograph, while FORS2 w as the instrumen t on the VL T. 
'DB' refers to objects from lDe Breuck et al.l (|2004 120061 ) . 

(8) Linear size in kpc using the largest angular size mea- 
surements at 2368 MHz unless marked with a footnote. 

(9) Rest-frame luminosity at 1.4 GHz in WHz -1 . Luminosi- 
ties were calculated using the 5-point spectral index from 
IPaperlj (table 2) unless marked by a footnote. 

2 NVSS J144932-385657 has one line with FWHM> 
3000 km s" 1 but it is faint and the error on the width is 
large; since the other lines are narrow we do not classify it as a 
QSO. 
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Of the 164 sou rces wit h K-band or SuperCOS- 
MOS counterparts in iPaper III spectra were obtained for 
52 objects listed in Table [2] including 46 new spec- 
tra and six that are common to th e SUMSS-NVSS 
sampl e and therefore have spectra in |Pe Breuck et al] 
(2006). NVSS J231727-352606 is in both the SUMSS- 
NVSS and MRCR-SUMS S samples, and h as a if-band 
magnitude limit listed in |Pe Breuck et al.l {2004) and a 
redshift from iDe Breuck et all |2006h ": our new, deeper 
limiting magnitude (see Paper is shown in the ta- 
ble. NVSS J103441-394957 and the combined source 
NVSS J023601-314204/J023605-314235 have very faint 
red identifications in SuperCOSMOS, and were retained in 
the sample because the identification was not located until 
we had high-resolution ATCA images. 

The spectroscopy has resulted in 36 confirmed redshifts, 
12 continuum spectra, one that remained undetected and 
three that are obscured by stars. The spectra are shown 
in Fig. [1] they were smoothed lightly for the plots, using 
(typically) a 3-pixel boxcar. 

Parameters of the emission lines were measured follow- 
ing pr ocedures similar to those described in lRottgering et all 
1)19971 ); these parameters are listed in Table [3] The line cen- 
tre and width (FWHM) were determined by fitting a Gaus- 
sian on top of a continuum; the line flux was obtained by 
direct summation of the flux above the continuum over a 
wavelength range four times the FWHM, except where such 
a width would also include neighbouring lines. No attempt 
has been made to de-blend lines. The location of the contin- 
uum was determined by eye separately for each line. This is 
the major source of error in the line ratios. The errors on 
the line fluxes and widths were estimated by comparing the 
measured values using different continuum levels; for some 
of the weaker lines, these errors are substantial. Because we 
extracted the spectra from only the central few rows of each 
object (see Section ^. 2. 4|l we have not in general included all 
the flux from each emission line, so the line ratios we derive 
are characteristic of the nuclear regions only. 

Thirty-three per cent of the redshifts obtained are in 
the range 2 < z < 3.5, corresponding to the epoch in which 
the star formation rate density of the Universe was at a 
maximum (see Section fT} . These objects will form part of 
our sample for follow-up observations to look at the cluster- 
ing properties and environments of massive galaxies at that 
epoch. 

Three of the possible identifications turned out to be 
M-stars. These are assumed to be misidentifications with 
foreground stars rather than radio sta r s, as t he surface den- 
sity of radio stars from lHelfand et"ahl (|l999T ) is much lower 
than our observed surface density of M-star radio identifica- 
tions. The fraction of contaminating stars is similar to the 
SUMSS-NVSS sample in which three out of 35 spectra were 
misidentified stars. The stars are not included in our analy- 
sis. Twelve sources (23 per cent) showed a continuous spec- 
trum with no emission lines and one was not detected, leav- 
ing 25 per cent of the galaxy spectra with undetermined red- 
shifts. This is a smaller fracti on than in the SUMSS -NVSS 
sample and the USS sample of lDe Breuck et all l|200lf) which 
found 28 per cent and 35 per cent respectively. Of our twelve 
continuum sources, one was a spectrum taken with the ANU 
2.3-m in which the S/N was insufficient to detect any lines 
and will be reobserved on a larger telescope. Six were NTT 



spectra which may need to be repeated on the VLT to im- 
prove the signal-to-noise ratio, and the remaining five were 
VLT spectra. The 'redshift desert' is the region 1 < z < 2.3 
where the strong [O ll]A3727 line has been shifted out of the 
spectral window, but Lya has not yet moved into the win- 
dow. The continuum sources may be in this redshift range 
or, alternatively, the Lya line may be self-absorbed and the 
other emission lines too weak to be detected. 



3.1 Notes on individual sources 

NVSS J02360 1-3142 04/J023605 314235 (combined 
source, see Pape rTlh : The spectrum shows no emission 
lines, but the redshift is secure based on stellar features. 
NVSS J103441-394957: The spectrum shows broad QSO 
emission lines and the radio source has a larg e linear s ize. 
NVSS J103615-321659: We proposed in IPaper ij that 
this may be a double source with very unequal lobes. While 
our higher resolution images did not show a double struc- 
ture as the fainter lobe was resolved out, our spectrum is 
that of a z = 2.136 galaxy. This high redshift supports our 
identification as it would be very unlikely for a z > 2 galaxy 
to coincidentally be as close as this to a USS radio source. 
NVSS J 105917 303658: The emission lines are broad, 
typical of a quasar and it has compact radio structure with 
a largest angular size of 1.0 arcsec. The C IVA1549 line shows 
strong self-absorption and the blue wing of the Lya line ap- 
pears to have been absorbed completely (see Fig. I14J1 . Ab- 
sorption at t he re dshift of the host galaxy was found by 

1 Baker et al.l ()2002h in almost all compact steep-spectrum 
radio-loud quasars. 

NVSS J 112641 -383950: The line at 7468A; (marked 
with an asterisk in Fig. [TJ is the A1305 O I+Si II line ap- 
pearing in second order. 

NVSS J140223-363539: The 8640-MHz image in 
iPaperlJ showed contours extending towards an object 

2 arcsec south-west of the 2368-MHz peak. A spectrum of 
this object has confirmed the redshift to be 2.796. While it 
does not align with the radio structure, it would be highly 
unlikely that this could be a separate high-redshift object 
that coincidentally is near the radio source. The NVSS cen- 
troid is 2.5-arcsec north of the claimed identification. That 
may indicate the presence of another faint lobe to the north 
which is resolved out in the ATCA observations, suggesting 
that it could have a wide-angled tail structure. For the time 
being we assume the high-redshift source to be the identifi- 
cation due to its redshift, but we cannot explain the offset 
from the expected position. 

NVSS J140854-382731: The redshift is based on a sin- 
gle line, which we identify as [Oll]A3727, and the marginal 
presence of other lines in the spectrum. 
NVSS J 141428-320637: Shows evidence of recent star 
formation, with H5 absorption characteristic of an A-star 
population. 

NVSS J144206-345115: A quasar, with broad lines and 
associated absorption (see Fig. I14p . 

NVSS J151020-352803: There is only one bright line in 
the spectrum, which is spatially extended over 3.6 arcsec, 
which we identify as Lya, based on the marginal presence 
of confirming lines. 

NVSS J151215-382220: No lines were identified, but the 
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Table 2. Spectroscopy results including redshifts and the derived linear sizes and rest-frame luminosities. Details of each 
column are given in Section 3. 
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a Two NVSS sources have been shown to be the components of one radio source. 

6 SuperCOSMOS UKR optical position. c 2MASS JSkrutskie et alj|2006h position. 

e 2MASS K magnitude in 8-arcsec aperture. 

n This source was observed by |Pe Breuck et al.l (2004) in the SUMSS-NVSS sample with a ii"-band non-detection. We 
have reobserved it and improved the limit on the ii'-band non-detection. 

1 These sources were also part of the SUMSS-NVSS sample. The ii'-band magnitudes and positions given are from 
iDe Breuck et al.l ll2004h . 

' This object is at the position where the redshift was measured , but it may be a for eground star (see Section [371} . 

m Spectral indices used in c alculating the lum inosities arc from Klamcr et al. (2006) and the linear sizes are derived from 

LAS measurements either in IDe Breuck et ahl d2004T) or iKlamer et al.l (120061 ') ■ 

n Linear size is based on the largest angular size at 4800 or 8640 MHz. Luminosity was calculated using the 7-point spectral 
index (including 4800- and 8640-MHz data) from iPaperU (table 3). 
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Figure 1. Spectra of sources with prominent features indicated. The source name, redshift and observing run are shown on top of each 
spectrum. The zero level is indicated by the horizontal dotted line. The positions of the vertical dotted lines indicate the predicted 
wavelengths of the lines at the redshift shown and not the wavelength of the fitted peak. Parentheses mark marginal detections of lines 
at the measured redshift and CR represents a cosmic ray. An asterisk shows the A1305 O I+Si II line appearing in second order. 
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Figure 1 - continued 
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Figure 1 - continued 
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Table 3. Emission line measurements. 



Source 




z 


Line 




■^obs 

(A) 




Adpwhm 
(km s J 


Flux 

(xlO erg s cm ) 


NVSS J103441- 


-394957 


1.832 ± 0.008 


C iv A1549 




4389.0 ± 


1.2 


6350 ± 700 


162 ± 40 








C in] A1909 




5368.4 ± 


9.7 


14900 ± 5000 


99 ± 25 


NVSS J103615- 


■321659 


2.136 ± 0.012 


Ly a A1216 




3803.0 ± 


0.3 


< 1700 


0.85 ± 0.21 








C iv A1549 




4871.3 ± 


0.5 


1580 ± 160 


0.30 ± 0.08 


NVSS J105917- 


-303658 


3.263 ± 0.022 


Ly 13 A1026 




4341.1 ± 


3.4 


7500 ± 1000 


1.6 ± 0.4 








Ly a A1216 




5219.8 ± 


1.0 


4200 ± 300 


1.9 ± 0.5 








C iv A1549 




6611.4 ± 


2.6 


7300 ± 600 


2.6 ± 0.6 








He II A1640 




6991.2 ± 


0.4 


1300 ± 80 


1.3 ± 0.3 








C in] A1909 




8124.8 ± 


0.9 


1600 ± 170 


1.5 ± 0.4 


NVSS J111921- 


-363139 


2.768 ± 0.006 


Ly a A1216 




4588.6 ± 


0.1 


270 ± 20 


1.8 ± 0.5 








C iv A1549 




5841.1 ± 


0.7 


2700 ± 180 


1.4 ± 0.3 








He II A1640 




6173.8 ± 


0.4 


1450 ± 90 


0.9 ± 0.2 








C in] A1909 




7179.3 ± 


0.6 


1200 ± 130 


0.9 ± 0.2 


NVSS J112641- 


383950 


1.856 ± 0.003 


C iv A1549 




4429.8 ± 


0.1 


750 ± 40 


0.6 ± 0.2 








He II A1640 




4687.4 ± 


0.2 


480 ± 60 


0.4 ± 0.1 








Al in] A1858 




5299.7 ± 


1.3 


1100 ± 400 


0.05 ± 0.03 








C in] A1909 




5448.4 ± 


0.1 


810 ± 30 


0.4 ± 0.1 








C II] A2326 




6647.2 ± 


0.4 


1270 ± 90 


0.4 ± 0.2 








[Neiv] A2423 




6921.2 ± 


0.5 


< 930 


0.15 ± 0.04 








Mg II A2798 




7998.0 ± 


0.3 


1120 ± 70 


0.9 ± 0.2 


NVSS J120839- 


-340307 


1.120 ± 0.002 


Mg n A2798 




5932.7 ± 


0.8 


1750 ± 450 


5.2 ± 3 








[0 n] A3727 




7906.7 ± 


0.4 


800 ±100 


7.5 ± 4 


NVSS J140022- 


-344411 


0.7844 ± 0.0007 


[O ll] A3727 




6652.4 ± 


0.3 


< 470 


7.6 ± 4 


NVSS J140223- 


-363539 


2.796 ± 0.012 


Ly a A1216 




4603.6 ± 


0.3 


1630 ± 100 


0.9 ± 0.2 








C iv A1549 




5892.1 ± 


0.3 


< 1100 


0.06 ± 0.03 


NVSS J140854- 


-382731 


0.870 ± 0.003 


[O n] A3727 




6966.4 ± 


0.2 


470 ± 90 


0.79 ± 0.4 








[Ne in] A3869 




7214.3 ± 


1.9 


1580 ± 700 


0.57 ± 0.3 


NVSS J141428- 


320637 


0.2775 ± 0.0008 


[O n] A3727 




4759.8 ± 


0.1 


660 ± 50 


125 ± 30 








[O in] A4363 




5576.2 ± 


0.1 


< 480 


5.7 ± 3 








H/3 A4861 




6209.9 ± 


0.1 


110 ± 50 


16.9 ± 4 








[O in] A4959 




6331.6 ± 


0.4 


560 ± 200 


13.2 ± 3 








[O in] A5007 




6393.5 ± 


0.1 


470 ± 50 


49.9 ± 10 








[N n] A6583 + 


Ha A6562 


8392.5 ± 


0.5 


2500 ± 100 


180 ± 45 


NVSS J144206- 


-345115 


2.750 ± 0.004 


Ly $ A1025 + 


O vi A1034 


3896.2 ± 


2.3 


8200 ± 900 


2.6 ± 0.6 








Ly a A1216 




4563.8 ± 


0.1 


2840 ± 50 


6.1 ± 1.5 








Si iv 1400 + O iv] A1402 


5252.5 ± 


2.3 


6100 ± 700 


1.1 ± 0.3 








C iv A1549 




5808.1 ± 


1.2 


6800 ± 300 


2.8 ± 0.7 








He ii A1640 




6153.0 ± 


9.2 


5800 ± 2000 


0.3 ± 0.1 








C in] A1909 




7142.3 ± 


1.0 


5600 ± 200 


1.8 ± 0.4 


NVSS J144932- 


-385657 


2.149 ± 0.002 


Ly a A1216 




3832.1 ± 


0.2 


2730 ± 70 


1.9 ± 0.5 








C iv A1549 




4876.0 ± 


0.2 


390 ± 50 


0.4 ± 0.1 








He ii A1640 




5164.6 ± 


6.4 


4700 ± 2000 


0.26 ± 0.06 


NVSS J151020- 


-352803 


2.931 ± 0.006 


Ly a A1216 




4782.2 ± 


0.1 


1340 ± 40 


1.4 ± 0.3 








C iv A1549 




6097.8 ± 


2.2 


1700 ± 500 


0.11 ± 0.03 








He ii A1640 




6436.9 ± 


0.9 


1330 ± 200 


0.11 ± 0.03 








C in] A1909 




7503.5 ± 


0.5 


1110 ± 100 


0.25 ± 0.06 


NVSS J152123- 


-375708 


1.530 ± 0.003 


C iv A1549 




3913.6 ± 


1.2 


8900 ± 400 


3.5 ± 0.9 








C in] A1909 




4831.5 ± 


0.6 


2600 ± 200 


0.6 ± 0.2 








Mg ii A2798 




7084.8 ± 


2.5 


5700 ± 500 


1.3 ± 0.3 


NVSS J152435- 


-352623 


2.752 ± 0.007 


Ly a A1216 




4571.5 ± 


0.3 


4430 ± 100 


3.1 ± 0.8 








C iv A1549 




5803.0 ± 


6.1 


15800 ± 1500 


2.5 ± 0.6 








C in] A1909 




7157.7 ± 


3.8 


7700 ± 800 


1.8 ± 0.4 


NVSS J152445- 


-373200 


1.798 ± 0.008 


C iv A1549 




4334.7 ± 


0.6 


1530 ± 400 


9.5 ± 5 








[Ne v] A3345 




9363.4 ± 


0.4 


740 ± 100 


20.8 ± 10 


NVSS J152737- 


-301459 


0.969 ± 0.008 


[O n] A3727 




7341.0 ± 


0.3 


340 ± 70 


0.61 ± 0.15 



brightness of the object and the presence of a significant very broad C IVA1549 line of 15800 kms" 1 FWHM confirm 
continuum make it unlikely to be at high redshift. that this is a QSO. 

NVSS J152123-375708: A QSO with broad ClvA1549. 

NVSS J213238-335318: The redshift is based on a single 
NVSS J152435 352623: Three broad lines including a emission line, which we identify as Lya, based on the absence 
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Adfwhm 


Flux 
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Line 


(A) 


(km s J 


/vyin — 16 —1 —2\ 

(xW ergs cm ) 


NVSS J210814- 


350823 


1.879 ± 0.002 


C IV A1549 


4462.1 ± 0.3 


1730 ± 100 


0.9 ± 0.2 








He II A1640 


4721.2 ± 0.2 


830 ± 70 


0.7 ± 0.2 








C in] A1909 


5491.1 ± 0.1 


< 1200 


0.5 ± 0.1 








C n] A2326 


6694.1 ± 0.6 


1340 ± 130 


0.24 ± 0.06 








[Ne iv] A2423 


6977.5 ± 0.4 


900 ± 80 


0.4 ± 0.1 








Mg II A2798 


8060.1 ± 0.4 


< 800 


0.4 ± 0.2 


NVSS J213238- 


335318 


2.900 ± 0.001 


Ly a A1216 


4740.6 ± 0.5 


2100 ± 200 


0.97 ± 0.2 


NVSS J213434- 


302522 


0.6819 ± 0.0005 


[O II] A3727 


6273.1 ± 0.2 


310 ± 90 


3.6 ± 2 


NVSS J213637- 


340318 


2.772 ± 0.002 


Ly a A1216 


4587.9 ± 0.6 


2700 ± 200 


1.0 ± 0.2 








C iv A1549 


5841.4 ± 1.3 


2850 ± 300 


0.44 ± 0.1 








He II A1640 


6184.1 ± 0.9 


2580 ± 200 


0.4 ± 0.1 








C in] A1909 


7201.6 ± 0.7 


1260 ± 140 


0.29 ± 0.07 


NVSS J215226- 


341606 


1.277 ± 0.002 


Mg II A2798 


6383.2 ± 0.4 


280 ± 100 


2.7 ± 1 








[Ne v] A3426 


7789.1 ± 0.7 


800 ± 200 


4.3 ± 2 








[O II] A3727 


8492.0 ± 0.9 


970 ± 300 


6.0 ± 3 


NVSS J215455- 


363006 


1.235 ± 0.001 


[Ne v] A3426 


7659.1 ±0.1 


< 410 


0.85 ± 0.4 








[O n] A3727 


8323.3 ± 0.4 


480 ± 100 


0.89 ± 0.4 


NVSS J223101- 


353227 


1.834 ± 0.006 


C iv A1549 


4400.7 ± 0.6 


1060 ± 200 


0.34 ± 0.08 








He II A1640 


4644.3 ± 1.5 


1300 ± 500 


0.14 ± 0.03 








C in] A1909 


5402.2 ± 0.8 


2060 ± 200 


0.6 ± 0.1 








C n] A2326 


6580.6 ± 1.4 


3690 ± 300 


1.0 ± 0.2 








[O n] A2470 


6995.8 ± 0.6 


1070 ± 120 


0.30 ± 0.07 








Mg n] A2798 


7945.5 ± 0.9 


2850 ± 200 


1.3 ± 0.3 


NVSS J223111- 


371459 


1.248 ± 0.002 


[O n] A3727 


8380.7 ± 0.4 


790 ± 100 


11.2 ± 3 


NVSS J223305- 


365658 


0.939 ± 0.002 


C II] A2326 


4519.2 ± 1.3 


1650 ± 800 


0.59 ± 0.3 








[Ne iv] A2423 


4691.9 ± 3.3 


2700 ± 2000 


0.67 ± 0.3 








[Ne v] A3426 


6643.5 ± 0.6 


670 ± 200 


0.28 ± 0.1 








[O n] A3727 


7222.2 ± 0.1 


680 ± 50 


1.39 ± 0.7 








[Ne in] A3869 


7500.9 ± 0.4 


520 ± 100 


0.40 ± 0.2 


NVSS J230004- 


304711 


0.5445 ± 0.0006 


[O n] A3727 


5755.3 ± 0.8 


1340 ± 400 


7.5 ± 4 


NVSS J232125- 


375829 


2.204 ± 0.006 


C iv A1549 


4972.8 ± 0.5 


2970 ± 300 


42 ± 10 








C in] A1909 


6103.4 ± 4.0 


5000 ± 2000 


11.5 ± 3 








Mg II A2798 


8975.8 ± 0.8 


830 ± 200 


5.1 ± 2 


NVSS J233034- 


330009 


2.675 ± 0.005 


Ly a A1216 


4461.7 ± 0.3 


960 ± 200 


4.4 ± 2 








C in] A1909 ? 


5687.8 ± 0.2 


< 550 


0.52 ± 0.3 








C iv A1549 ? 


7017.2 ± 3.8 


4200 ± 1500 


3.9 ± 2 


NVSS J233226- 


363423 


0.988 ± 0.003 


Mg n A2798 


5576.1 ± 10 


10700 ± 500 


19.2 ± 5 








[O n] A3727 


7400.1 ± 1.3 


2000 ± 500 


11.2 ± 3 
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Figure 1 - continued 
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of blueward continuum emission. The bright line is spatially 
extended over 1.8arcsec, supporting the Lya identification. 
NVSS J223101-353227: The Cll]A2326 and MgllA2798 
lines are broad, so we classify this object as a QSO. The line 
ratios support this interpretation (see section I4.5|l , though 
the line equivalent widths are unusually small for a QSO. 
NVSS J223111-371459: The spectrum shows both 
[Oll]A3727 emission and strong Balmer absorption indica- 
tive of an A-star population. 

NVSS J232125 375829: A quasar. The blue end of the 
spectrum is noisy and has been truncated but Lya is not 
detected. 

NVSS J233034 330009: From our original 1384- and 
2368-MHz radio images, we initially chose the identification 
to be the bright Tf-band object between the lobes and closer 
to the stronger lobe. We had taken a spectrum of this ob- 
ject befo re the 4800-MHz-high-resolution image (shown in 
IPaper ij) identified a possible core which is 4-arcsec east of 
the bright K-band object. The spectrum is typical of an M- 
star, but there is a clear emission line in the blue which we 
identify as Lya at a redshift of 2.675. The slit position an- 
gle was 165° and centred on the bright object and therefore 
does not go through the core position. The spatial profile of 
the Lya line, shown in Fig. [2j is extended with a FWHM of 
7.4arcsec (60kpc) along the slit. We therefore believe that 
the correct K-band identification is fainter than the limit 
on the IRIS2 image (i.e. K > 19) and that it is surrounded 
by a Lya halo which extends from the core to beyond the 
bright source, which is an M star. In Fig. [3] we show the 
possible position of the Lya halo. In a simplistic picture, 
where the Lya halo is spherical, the 7.4-arcsec FWHM of 
the emission line at the slit can be reproduced by a halo of 
diameter ~ 10.5 arcsec, or 85kpc. With an alternate geom- 
etry, the halo size would need to be a minimum of 32kpc, 
the projected distance from the core to the bright object. 
This is a modest size f or a Lya halo, as some have been 
found to be > 100 kpc jKurk et al.ll2000al ; IVenemans et all 
120021 ; iReuland et alj|2007h . Since the possible core was not 
detected in our 8640-MHz image this interpretation is open 
to question. 

An alternative picture is that the bright K-band object 
is a galaxy and star coincident along the line of sight, and 
the middle component is a hotspot along the eastern jet. 
If this picture is correct then there is a marginal line at 
the wavelength expected for C IVA1549. A further line sits 
at the expected position of Cm]A1909, but it cannot be 
differentiated from the M-star features at that wavelength. 
NVSS J233226 363423: A very broad MgllA2798 line, 
with a FWHM of 10 4 kms" 1 indicative of a QSO. 




Offset (arcsec) 

Figure 2. Spatial profile (north-south, positive offsets to 
the north) of the Lya line in the NTT spectrum of 
NVSS J233034- 330009. Lya covers a FWHM of 7.4 arcsec or 
60 kpc. 
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4 DISCUSSION 
4.1 K—z distribution 

4-1.1 Emission-line contribution to K-band magnitudes 

Strong emission lines such as [O 11] A3727 , [O in] AA4959, 5007 
and Ha A6562 can be shifted into the A"-band, particularly 
at redshifts beyond ~ 2. The contribution of these lines 
to the if-band magnitude is primarily dependent on red- 
shift, but also on radio luminosity, due to t he correlation 
between emission-line and radio luminosities (jWillott et alj 



Figure 3. Our proposed model for NVSS J233034-330009 (see 
Section [311 drawn on the overlay of the 4800-MHz ATCA con- 
tours on the smoothed ii"-band image. The lowest contour is 3 
sigma, and the contours are a geometric progression in \/2. The 
rms noise (er) is shown in the header in mjybeam -1 . The par- 
allel lines represent the slit orientation (position angle of 165°) 
through the bright source (with cross- hairs), assumed to be a star. 
The dashed circle is the proposed Lya halo, centred on the radio 
core, and with a diameter chosen to reproduce the observed width 
of the Lya emission at the slit (7.4 arcsec), giving a halo diameter 
of 10.5 arcsec or 85 kpc. 
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1 19991 ). |j arvis et all l|200lf ) calculated the emission-line con- 
tamination at each redshift for different 151-MHz flux den- 
sities. Adopting their approach, we have estimated the 
emission-line contamination to our A-band magnitudes for 
the MRCR-SUMSS sources. Firstly, we extrapolated our ra- 
dio spectrum to 151 MHz, using the spectral index from our 
5-point fit or, whe n unavail able, using the 843-1400-MHz 
spectral index (see IPaper III) . Some of the sources have a 
measured 151-MHz flux density (see Table [5]) which was 
used. We then estimated the A-band contribution from 
emission lines using the information from ljarvis et all (|200ll . 
fig. 2), and then subtracted the emission-line contribution 
from the measured A flux. We find that in all but one 
case (NVSS J144932-385657) the typical contribution of 
the emission lines is less than the errors on our K magni- 
tudes. The inherent assumptions in this process mean that 
a typical emission-line contamination was calculated, rather 
than an actual emission-line contribution. These assump- 
tions include using average line ratios, the generalisation of 
the correlation between emission-line and radio luminosity, 
and the assumption that our radio spectra remain straight 
down to 151 MHz. 

The only source where line emission may contribute to 
the K flux is NVSS J144932-385657 at z = 2.149, where 
the emission lines of [N n]A6583 and Ha fall in the K band. 
This source is relatively faint at K — 19.8 ± 0.2 and the 
possible emission-line contribution may account for up to 
0.4 magnitudes of A-band flux. However, the spectrum in 
Fig. [1] shows that the emission lines are weak and therefore 
the lines that are shifted into the A-band may not have the 
maximum flux calculated here. Furthermore, this object is 
not an outlier in the K-z distribution (see Fig. [4] below) so 
there is no evidence that emission lines have boosted its K 
magnitude. 



4-1.2 Trends on the K-z correlation 

The K-z correlation is an empirical relationship that spans 
redshifts up to z > 4. While the correlation remains lineai0, 
different galaxy samples have b een shown to give d ifferent 
linear fits to the K-z relation. IWillott et ail l|2003h made 
a fit to 64-kpc-aperture magnitudes based on the 3CRR, 
6CE, 7C and 7CRS galaxy samples. The 3CRR sample was 
found to be ~ 0.6 magnitudes brighter in K at z > 1 
than the 6C l|Jarvis et all 120011 ). while the 7 C and 3CRR 
are o ffset in A-magnitude at all redshifts (|Willott et al.l 
120031). A fit to the K-z relation for FIRST radio sources 



1 E1 Bouchefrv fc Cress 2007) was also found to be fainter 
than the Willott et al.l (120031 ) relation by more than 0.5 mag- 
nitudes at high redshift and even more at low redshift. These 
differences highlight the selection effects (including flux den- 
sity limit) in any sample and show that the fits to the K-z 
plot depend on which samples are included and the aperture 
correction used. An individual fit, therefore, has no physical 
meaning but serves only to highlight whether a galaxy lies 
above or below the trend for that sample. 

In Fig. U the 4- arcsec- aperture K magnitudes are plot- 
ted against redshift for all the sources in our sample that 

3 Linear fit between log(z) and K magnitude corresponds to a 
power-law relation between z and 2.2fim flux. 



have spectroscopic redshifts, along with many samples from 
the literature with 4-arcsec-aperture K magnitudes. These 
are both USS-selected samples and non-USS-selected sam- 
ples as liste d in the Figure, and include the largest non- 
USS sample l|McCarthv et al.lll996l ; 2007, private communi- 
cation) comprising 277 A-band observations of radio galax- 
ies, of which 175 have redshifts. The advantages of using 
the 4-arcsec-aperture magnitudes are discussed in IPaper III 
and in terms of the K-z dispersion in Section ^. 1.31 We have 
therefore chosen to adopt a fit to the 4-arcsec-aperture mag- 
nitudes from the combined literature samples as the most 
suitable metric that represents our data and allows more of 
the galaxies in our sample to be compared with the linear 
relation for HzRGs. For consistency with some papers in 
the literature we also show the K-z plot based on 64-kpc 
aperture values. The MRCR-SUMSS 64-kpc magnitudes are 
calculated from 8-arcse c-aperture magnit udes following the 
procedure described bv lEales et alj (|l997l ) and listed in Ta- 
bled 

Three fits to the 4-arcsec-aperture magnitudes in Fig- 
ureUare shown. The fits are linear least-squares fits, because 
in each case the error on a quadratic fit was higher. QSOs 
and objects with contamination from foreground stars were 
not included in the fits. The equations of the fit lines are 

A(4") = 3.64 log(z) + 17.75 (1) 

for objects at all redshifts (solid line), 

K(4") = 3.45 log(z) + 17.76 (2) 

for z > 0.6 objects only (dotted line), and 

A(4") = 3.11 log(z) + 17.89 (3) 

for z > 1.0 objects (dashed line). 

As the redshift range increases to include lower redshift 
objects, the fit steepens, whereas if the 4-arcsec aperture 
was missing any A'-band flux we would expect the opposite 
effect, i.e., the fit would be flatter. There is therefore no 
compelling evidence that the lower redshift galaxies have 
A-band emission detected beyond a 4-arcsec aperture. 

In Fig. 2] we have also plotted the evolutionary 
curves for elliptical galaxi es with masses of W^Mp an d 
1O 12 M from the models of lRocca-Volmerange et all l|2004h . 
The majority of identifications in our sample are con- 
sistent with elliptica l galaxies of mass 10 11 - 10 12 M . 
ISevmour et all i|2007h modelled HzRG masses from H— band 
magnitudes using early-type gala xy models from PEGASE2 
|Fioc fc Rocca-Volmerangall997l ) with galaxy formation at 
z = 10 and a different init i al mas s function to th a t used 
by iRocca-Volmerange etail (|2004l ). ISevmour et all (|2007f) 
found that the majority of their galaxies fitted within the 
10 11 -10 11 ' 5 M Q model over a redshift range of 1 < z < 4. 

The K-z distribution has a scatter about the fitted line 
due to the variations in stellar content of the galaxies, the 
amount of dust, and the rest wavelength range being sam- 
pled at each redshift. At high redshift UV light is observed 
in the near-IR, and this may introduce a contribution from 
the AGN continuum. In Fig. [3] there are two galaxies that 
are more than 2.5<r from the adopted best fit line (Eq. [IJ, 
and lie well outside the 10 11 Mq and 10 12 Mq model lines. 
NVSS J011606-331241 is 3.3cr fainter than the fit to the 
main distribution of points, with K — 18.6 and z = 0.352. 
This galaxy is in common with the SUMSS-NVSS sample; 
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Figure 4. Top: K magnitudes in 4-arcsec apertures versus redshift for the MRCR-SUMSS HzRGs (filled circles), the MRCR-SUMSS 
QSOs (stars) and samples from the literature (small open circles). The literature samples were chosen because they had 4-arcsec 
iC-band and redshift information. QS Os, uncertain redshif t s, an d lower limi t K magnitudes were removed from the literature sam- 
ples. Two of the samples (CENSORS. iBrookes et al.ll2006l . 120081 : and 6C**. ICruz et aljliopj l had 3- and 5-arcsec-aperture K mag- 
nitudes and a 4-arcsec-aperture magnitude was calculated by av eraging the equival e nt 3- and 5-arcsec fluxe s and converting to a 4- 
arcsec magnitude. T he other samples plo tted are the 6C* from Uarvis et al.l (200 it ), foe Breuck et al 1 ll200ll. |2002i) sample, 7C from 



IWillott et ahl 11200311. SUMSS-NVSS from foe Breuck et al.l d200^ . l200d) andthe sample from iBornancini et al.l ||2007|) . The small dots 
are the " [McCarthy et al. (1996; 2007. private cmr!rminication)l sample in which the magnitudes were measured in a 3-arcsec aperture. 
Three straight line fits to both the literature samples and the MRCR-SUMSS sample are shown; a fit to points at all redshifts (solid 
line, Eq. [TJ, to points with 2 > 0.6 (dotted line, Eq. [2j and to points with z > 1 (dashed line, Eq. [3}. The McCarthy points were 
not used in the fit calculations due to the smaller measurement aperture for that sample and neither were the QSOs and objects with 
contamination by a foreground star in the MRCR-SUMSS sample. The redshift errors for our MRCR-SUMSS data are much smaller 
than the plotted filled circle as are some of the X-magnitude errors. Bottom: K magnitudes in a 64-kpc metric aperture versus redshift 
for the MRCR-SUMSS HzRGs (filled circles), the MRCR-SUMSS QSOs (stars) and sample s from the lite r ature (small open circles). The 
lit erature samples ar e the 3 CRR , 6CE, 6C* and 7CRS radio galaxy samples compiled bv lWillott et al ] l|2003t ). the composite samples 
of Ivan Breueel et al.l lll998ll and foe Breuck et al.l ||2002|) , the SUMSS-NVSS sample and the McCarthy sample. The 64-kpc aperture 
magnitudes for the MRCR-SUMSS, SUMSS-NVSS and McCarthy samples were calculated from 8-arcsec- aperture magn i tudes (see text 
for limit ations). We have a lso pl otted the corresponding 64-kpc and 8-arcsec-aperture fit equations from lwillott et al.l J2003T) (dashed 
line) and |DeBreuck_e^^lJ (2004) (dot ted lin e) respectively. The dot-dash lines are the 10 11 Mq and 10 12 Mq elliptical galaxy evolution 
models from lRocca-Volmerange et al.l 1(2004]). 
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Table 4. The dispersion in the X-magnitude as a function of 
redshift for all the galaxies in Fig. [4] but excluding QSOs and 
galaxies that are contaminated by foreground stars. The number 
of galaxies in each redshift bin is given in brackets beside the 
dispersion. 



Fit equation 


Dispersion (mag) 




all z's 


z < 2 


z > 2 


64-kpc fit 


0.7 (534) 


0.7 (441) 


0.9 (93) 


all 4" ap. 


0.7 (211) 


0.7 (146) 


0.7 (65) 


z > 0.6 for 4" ap. 


0.7 (211) 


0.8 (146) 


0.7 (65) 


z > 1.0 for 4" ap. 


0.8 (211) 


0.8 (146) 


0.7 (65) 



in iDe Breuck et all (|2006h it was concluded that there is a 
chance the correct identification may be a fainter nearby 
galaxy but higher resolution radio imaging is needed to con- 
firm this. The other main outlier is NVSS J233034-330009 
which is 2.8cr brighter than predicted, with z — 2.675 and 
K = 17.2; this source is discussed in detail in Section T3, II 

4-1.3 Dispersion in the K-z distribution 

An increase in the dispersion of points on the K-z plot 
at high redshift has been used in the past to co n strain 
the epoch of formation of radio galaxies. lEales et al.l (|l997l ) 
found an increase in dispersion above z ~ 2, and concluded 
that those galaxies were in the process of formation while 
those at z < 2 ar e evolving passively. On the other hand, 
Ijarvis et alj (|200ll ) found no evidence for an increase in the 
dispersion of K magnitudes at z > 2, and concluded that 
radio-luminous systems formed most of their stars at z > 2.5 
and have evolved passively since then. 

In Table [4] we have calculated the dispersion in redshift 
bins for all the galaxies in Fig. [3] based on each fit shown. 
This shows that the disper s ion in the 64-kpc-aperture points 
around the IWillott et alj (|2003T > fit is no lower than dis- 
persion of the 4- arcsec- aperture points about the fit to all 
redshifts. Therefore, we have adopted the 4-arcsec-aperture 
metric with the fit in Eq.[T]in order to compare the MRCR- 
SUMSS galaxies with the trend in radio galaxies for the 
rest of this paper. There is no evidence for an increase 
in dispersion above z = 2 for th e 4-arcsec-aperture fits. 
This supports the resu lts found bv lJarvis et alj (|200ll ) and 
Dc Breuck et al.1 (|2006j ) that radio galaxy hosts have been 
evolving passively since epochs corresponding to z = 3 or 
earlier. 

4.1.4 QSOs in the MRCR-SUMSS sample 

QSOs are usually removed from HzRG samples on the ba- 
sis that their if-band flux is generated by the central AGN 
rather than from stellar light. The implication is that the 
non-thermal AGN flux would be stronger than the ther- 
mal contribution from the host and hence bias the K— 
z distribution in favour of brig hter if -band magnitudes. 
The unification model of AG Ns l|Barthell [19891 ; lAntonuccil 
1 19931 : lUrrv fc Padovanil Il995l ) suggests that the difference 
between radio galaxies and radio-loud QSOs is an orien- 
tation effect, such that the QSOs are viewed so the nu- 
cleus can be seen directly rather than being obscured by 



a dusty torus as in a radio galaxy. QSOs would therefore 
be expe c ted t o have smaller apparent linear sizes. While 
iBarthell 0989) found that QSOs, on average, were a fac- 
tor of two smaller than radio galaxies, many QSOs have 
large radio linear sizes, up to several Mpc l|Rilev et all 
1 19891 : iBhatnagar. Krishna fc Wisotzkil 11998 : iPaperlll se~ 
tion 3.2.1) and some do not show very bright central sources. 
Alternatively, it has been argued that the difference between 
QSOs and radio galaxies may be due instead to their evo- 
lutionary stage and/or environment. For instance, compact 
radio sources were fou nd to live in denser environm ents than 
extended sources by iPrestage fc Peacock! (|l988l ) and this 
was used as an argument for them being different types of 
sources. It is possible, therefore, that radio morphology may 
depend more on the age or environment of a source than its 
viewing angle. As a result we consider it important to in- 
clude QSOs in any follow-up study of HzRG environments. 

It has been suggested that the alignment effect be- 
tween radio structures and UV /optical light se en in HzRGs 
is also apparent in QSOs (|Lehnert et al.lll999l ) and can be 
explained by scattering of light from the active nucleus 
rather than jet-induced star formation. Scattering models 
are favoured because the emission is highly polarised. While 
broad emission lines have long been the defining features 
of QSOs, the distinction between radio-loud QSOs and ra- 
dio galaxies is not clear-cut. It is worth noting that com- 
pact steepj^p_e£tnrm_QSOs in the Molonglo Quasar Sam- 
ple (|Baker fc Hunste ad 1995) also have properties similar to 
HzRGs: a heavily reddened optical continuum, strong nar- 
row forbidden lines and narrower than usual permitted lines. 
USS QSOs are quite rare, and therefore may have different 
physical characteristics to the bulk of the QSO population. 

While our fits to the K-z plot do not include the QSOs, 
it can be seen in Figure[4]that the QSOs tend to fall on the 
faint side of the fit to all redshifts, indicating that the K- 
band flux is not unusually strong in comparison with the 
radio galaxies. If the QSOs had been included, the 2 > 1 fit 
does not change Eq. [3]within the precision shown. We note, 
however, that target selection for spectroscopic followup was 
based on 7f-band magnitude but not on Jf-band morphol- 
ogy; hence bright QSOs will not be included in our sample. 

While the objects with line widths > 3000 km s~ have 
been classified as QSOs, we find no evidence from the K- 
band magnitudes that we are seeing the nuclear source di- 
rectly. Furthermore, two of the MRCR-SUMSS QSOs - 
NVSS J223101-353227 and NVSS J233226-363423 — have 
extended if-band emission that is aligned (within 7- and 3- 
degrees respectively) with the radio structure. 

In Fig. [S] we compare the MRCR-SUMSS radio galax- 
ies and QSOs with the overall K-z fit for radio galaxies de- 
fined by 4-arcsec aperture magnitudes. The offsets between 
the measured redshifts and those predicted by Equation [T] 
z(fit)— z(spec), are plotted against radio luminosity, linear 
size and redshift. It is clear that the differences among the 
K-z fits defined by Eqs. [T] [5] and [3] are small in compari- 
son to the scatter about each fit. Moreover, these QSOs do 
not stand out from the HzRGs in linear size, radio luminos- 
ity or redshift distribution. It is only the spectral line width 
that differentiates QSOs from HzRGs in the MRCR-SUMSS 
sample. 
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Figure 5. The difference between predicted and spectroscopic rcdshifts (where the predicted redshift comes from the K—z fit in Fig. [4j 
Eq[TJ as a function of linear size, 1.4-GHz luminosity and spectroscopic redshift. Circles are radio galaxies and crosses mark the QSOs. 



4.2 Linear size 

In some USS-selected HzRG samples an explicit maximum 
rad io angular size cut-off has been applied (see f or exam- 
ple iBlundell etlll Il998l ; |Pe Breuck et all l2000ah . on the 
assumption that the higher-redshift sources are yo unger 
and hence smaller (IBlundell. Rawlings fc Wi llott 1999). Our 
plot of largest linear size versus redshift in Fig. HI] is limited 
by the small number of redshifts measured so far (21 per cent 
of K detections have redshifts) . The selection of galaxies for 
redshift followup is likely to exclude some radio-extended 
but optically-faint objects (as discussed in Section I4.4.2p 
which we expect to be lower redshift dust-obscured galax- 
ies that would populate the top left of Fig. [5] However, 
the upper envelope of source sizes drops at higher red- 
shift. Linear size evolution is well known in radio samples 
dKapahi. S ubrahmanva & Kulkarnl Il987l ; iBarthel fc Milevl 
1 19881 ), but IBlundell et all (|2002h found only a weak cor- 
relation between redshift and linear size measured at low 
(74 MHz) frequency. The correlation here is also weak. Fur- 
ther analysis will be carried out when the spectroscopy for 
the MRCR-SUMSS sample is complete. 

The measured if —magnitudes have also been plotted 
against largest angular size (LAS) in Fig. [7] Our adopted 
cosmology dictates that the LAS will change with redshift, 
and therefore a line is shown representing the LAS of a 
50-kpc source at each redshift, and similarly for a 100-kpc 
source. In both cases, the LAS was calculated using our cho- 
sen cosmological parameters (see Section [1} , and then the 
redshifts were converted to K magnitudes using Equation [T] 
If the distribution of points decreases faster than the lines 
marking the cosmological models, then more distant USS- 
selected galaxies may be intrinsically smaller than their low- 
redshift counterparts. The median LAS of our objects at 
each K magnitude generally decreases faster than the cos- 
mological models. However, a Spearman rank test confirms 
there is no correlation in Fig. [7] 



4.3 The alignment effect 

Alignment between near-infrared and radio structures has 
been found in many H zRG samples, and these are discussed 
in detail in iPaper I J . However, the redshift dependence of 
the alignment effect is very much in dispute, with some 
investigation s finding a clear increa se in alignments with 

while others find no 
We now examine the 



redshift (e.g. Ivan Breugel et al" 
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6. Largest linear size from Table [2] versus redshift for 
(excluding QSOs) in our sample with spectroscopic red- 




K— magnitude (4 arcsec ap.) 

Figure 7. Largest angular size (LAS) versus 4-arcsec-aperture K 
magnitudes for the MRCR-SUMSS sample. T he LAS w as mea- 
sured from the 4800- or 8640-MHz images (see lPaperlJ. table 3) 
where available, and otherwise, the value comes from IPaper I j . 
table 2. Two solid curves represent the angular size of a 50-kpc- 
and 100-kpc-diameter galaxy versus redshift (based on Equation 1 
from the K—z plot) for our adopted cosmological model. The heav- 
ier solid line joins points which are the median LAS values in each 
magnitude bin from 16—17 up to 20-21 magnitudes. Solid points 
represent the objects that were followed up with spectroscopy (ir- 
respective of whether that resulted in a redshift measurement). 
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Figure 8. The offset between the radio axis position angle and 
that of the major axis of the ellipsoidal fit to the A"-band galaxy 
versus redshift. Sources with LAS<5 arcsec are not included as 
the position angle offsets may be less accurate. 



redshift dependenc e of the a lignment effect in the MRCR- 
SUMSS sample. In lPaperlJ we fitted ellipsoidal profiles to 
the _K"-band host galaxies. The major axis of the ellipsoids 
was found to align with the radio structure in 21 per cent 
of the 122 galaxies for which an elliptical profile could be 
fitted. Fifteen of the 122 fitted galaxies have redshifts and 
they are plotted in Fig. [8] against the position angle off- 
set, defined as the difference between the radio axis and the 
major axis of the fitted 7\"-band ellipsoid. There may be an 
upper envelope to the distribution of points, which falls with 
increasing redshift. If that remains the case when the spec- 
troscopic sample is complete, then it would imply that the 
higher- redshift objects show closer alignment than those at 
lower r edshifts. 

In lPaperHl we identified sources with alignment offsets 
of less than 10 degrees as being likely candidates for jet- 
induced star formation. Only three of those have redshifts 
in Fig. [8] While the trend with redshift identified here is 
clearly not associated with such highly aligned processes, 
larger position angle offsets have been attributed to ionisa- 
tion cones or extended emission-line regions in some high 
redshift galaxies. Such features would not be seen in galax- 
ies at 1.3 ^ z ^ 1.9 where there are no bright emission lines 
in the Jf-band. It is noteworthy that the objects which are 
highly misaligned, including those with perpendicular axes 
(> 80 degrees), are all at low redshift. 



4.4 The efficiency of low-frequency selection 

4-4-1 Spectral curvature 

The original purpose of selecting sources that are ultra- 
steep between 408-843 MHz was to test whether the sam- 
ple would preferentially select higher redshifts than USS 
samples selected at higher frequencies. This would be the 
case if the radio spectrum steepens at higher frequencies 
and the steeper part is redshifted down to 408-843 MHz, 
in which case the lower-freq uency se lection would include 
fewer low-redshift objects. In IPaper il we found that 85 per 
cent of USS-selected sources have straight radio spectra be- 
tween 408-2368 MHz (observed frame) and therefore the fre- 



quency of USS selection should not affect the resultant sam- 
ple if selected at a fr equency higher than 408 MHz. Similarly 
Kl amer et al . (2006) found that 89 per cent of a subsam- 
ple of 37 sources from the SUMSS-NVSS sample also have 
straight radio spe ctra from 8 43 MHz to e ither 6.2, 8.6 or 
18 GHz. However, H ornancmi et al l l|2007h found that the 
majority of their USS sample have radio spectra that flat- 
ten from 352 MHz to 74 MHz. If the radio spectrum remains 
straight over a larger frequency range, then the spectral in- 
dex is independent of the observed frequency. In that case, 
the z—a correlation would not be due to a k-correction. The 
steep spectrum would instead be intrinsic to the source or 
the environment it is in. 

We have an opportunity to investigate the shape of the 
radio spectrum below 408 MHz for the MRCR-SUMSS sam- 
ple as some of the sources are al so in the 74-MHz V LA Low- 
Frequency Sky Survey (VLSS; ICohen et al.l [20071 ) and the 
151-MH z Mauritius R adio Telescope (MRT) Southern Sky 
Survey (|Pandev Il200?3 ). These two catalogues, while useful, 
have limitations. Firstly, the declination of our sources is 
at the southern edge of the VLSS survey region, where the 
flux densities are considered less reliable. Secondly, the large 
MRT beamwidth (4x4.8arcmin at declination —35°) is sub- 
ject to beam confusion. We note that the MRT flux density 
calibration is tied to i nterpolated 151-MHz flux de nsities for 
strong sources in the iBurgess fc Hunsteadl ||2006T ) compila- 
tion. Errors in flux density were assumed to be the quadratic 
combination of a 0.26-Jy-rms no ise plus confusi on term and 
a 6.3 per cent scale uncertainty i|Pandevll20"06r i. 

Inclusion of 74- and 151-MHz flux densities in our spec- 
trum fits was done only after careful assessment of the ra- 
dio images. MRT and VLSS images were overlaid with im- 
ages from MRCR and SUMSS to check for confusion in the 
MRT images, and flux density or morphological anomalies in 
VLSS at these large zenith distances. The MRCR-SUMSS 
sample sources which have a reliable 74-MHz flux density 
have been fitted with a 6 -point spectral index between 74 
and 2368 MHz in IPaper J (table 10) and those with a K- 
band counterpart are listed in Table [6] The 151-MHz cata- 
logue does not cover all of our survey region, but we found 
matches to ten of the MRCR-SUMSS sources which have K- 
band counterparts. These are given in Table [5] along with a 
revised spectrum fit which includes the new frequency point. 
Two of these have both 74-MHz and 151-MHz flux densities. 
An example of the fit to one source is shown in Fig. [9] 

In Tables [5] and [6] 75 per cent (15/20) of sources have 
spectral shapes that remain a power-law down to either 74 or 
151 MHz. The mean if-band magnitudes are 17.9 and 18.0 
for the flattening and straight-spectrum sources respectively. 
The difference between these distributions is not statistically 
significant (based on a Wilcoxon rank-sum test). Therefore, 
we find that the majority of our sources have a spectral 
shape that remains straight, which supports a picture in 
which the steep spectrum is intrinsic to high redshift sources 
and that a k-correction is not the reason why these high- 
redshift sources have been identified by our selection criteria. 



4-4-2 Relationships among redshift, K, a and 1400MHz 
flux density distributions 

In lPaperlll we compared the A'-band magnitude distribution 
of the MRCR-SUMSS sample to both USS-selected- and 
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Table 5. Sources from the MRCR-SUMSS sample which have a _R"-band identification and a 151-MHz flux density in the MKT catalogue. 
K magnitudes in a 4-arcsec aperture and the 151-MHz flux densities are listed. The spectral index is shown for sources with a power-law 
radio spectrum. All the spectral indices are 6-point values between 151 and 2368 MHz, except for NVSS J213238— 335318 which has an 8- 
point index from 151 to 8640 MHz, NVSS J233729-355529 which has a 4-point index from 151 to 1400 MHz and NVSS J234145-350624 
which has fluxes measured at 9 frequencies from 151 MH z to 18 G Hz. Sources without spectral indices listed flatten towards lower 
frequencies. Objects which also have 74-MHz flux densities (jPaper j , table 10) are marked with an asterisk, and the 74-MHz flux density 
was used in the spectral index fit. 



Source 




K 


Sl51 


a 


Redshift 






(4-arcsec) 


(Jy) 






NVSS J003445- 


372348 


17.0 


2.77 ± 0.31 


-1.20 ±0.03 




NVSS J202720- 


341150 


> 20.3 


2.39 ±0.30 


-1.13 ± 0.04 




NVSS J213238- 


-335318 


19.8 


1.32 ±0.27 


-1.41 ± 0.02 


2.900 


NVSS J214114- 


332307* 


19.1 


2.19 ±0.29 


flattening 




NVSS J215455- 


363006 


17.7 


3.78 ±0.35 


-1.13 ± 0.03 


1.235 


NVSS J223305- 


365658 


16.6 


1.99 ±0.29 


-0.95 ± 0.04 


0.939 


NVSS J233535- 


343330 


18.9 


3.40 ±0.34 


-1.05 ± 0.03 




NVSS J233729- 


355529 


19.2 


2.58 ±0.31 


-1.39 ± 0.05 




NVSS J234145- 


350624 


16.9 


11.26 ± 0.76 


flattening 


0.644 


NVSS J235945- 


330354* 


19.1 


2.76 ±0.31 


flattening 





Table 6. Sources from the MRCR-SUMSS sample which have a 74-MHz flux density f from lPaper J . table 10) and a A"-band identification 
and are not already shown in Table [5] K magnitudes in a 4-arcsec aperture and the 74-MHz flux densities are listed. All the spectral 
indices are six-point values between 74 and 2368 MHz except for NVSS J230846— 334810 which has a four-point fit between 74 and 
1400 MHz. Sources without spectral indices listed flatten towards lower frequencies. 



Source 


K 


S74 


a 


Redshift 




(4-arcsec) 


(Jy) 







NVSS J000231- 


342614 


18.2 


1.35 ±0.19 


-0.98 ± 0.04 




NVSS J000742- 


304325 


18.9 


0.77 ± 0.10 


flattening 




NVSS J001506- 


330155 


17.6 


2.10 ±0.24 


-0.89 ± 0.03 




NVSS J021759- 


301512 


18.3 


1.56 ±0.17 


-1.04 ± 0.03 




NVSS J120839- 


340307 


18.0 


22.93 ± 2.34 


-1.10 ± 0.03 


1.120 


NVSS J141428- 


320637 


15.5" 


3.78 ± 0.49 


flattening 


0.2775 


NVSS J215047- 


343616 


17.4 


4.49 ± 0.53 


-1.05 ± 0.03 




NVSS J215226- 


341606 


18.3 


1.67 ± 0.22 


-0.99 ± 0.03 


1.277 


NVSS J221650- 


341008 


16.4 


2.23 ±0.29 


-1.00 ± 0.03 




NVSS J230846- 


334810 


17.0 


1.61 ± 0.20 


-1.13 ±0.04 





a 2MASS magnitude in an 8-arcsec aperture. 



non-USS-selected-samples from the literature, and found 
that the resulting i^-band distribution was independent of 
the flux density distribution. Here we look at the redshift 
distribution of the same samples and we investigate any bi- 
ases in the selection of objects for spectroscopic followup. 

Fig. [TD] shows the redshift and 1400 MHz 
flux density distributions for our MRCR- 
SUMSS! sample an d the SUMS S- NVSS , 6C* * 
dCruz et all I2006TI, iDe Breuck et~aLl (l200ll. [2 002). 



iMcCarthv et al. (1996; 2007. private communication)! . 
CENSORS l|Brookes et all l2006l . |200S|) and a combined 
sample from the 3CRR, 6C*, 6CE and 7CRS catalogues 
(compiled bv lWillott et~al]|2003l ) . The McCarthy and CEN- 
SORS samples were not USS-selected, while the Willott et 
al. sample is predominantly non-USS-selected apart from 
23 sources from the 6C*. The McCarthy sample includes 
175 redshifts from a set of 277 if-band observations of 
radio galaxies. It is the largest non-USS K-band sample 
even though it remains unpublished. The 1400 MHz flux 
density values for the McCarthy sample were obtained 



by cross-matching the positions with the NVSS catalogue 
within a 60-arcsec search radius. For sources which had 
two matches within 60 arcsec, we selected the closest 
match. SUMSS-NVSS was selected to have < -1.3, 

and the 6C** sample is from the 151-MHz 6C surve 
with aif? < -1.0. The IDe Breuck et alj l|200ll . 1200 



sample was selected at one of the frequencies 325, 365, or 
408 MHz as detailed in the papers, and h as odj§Sn S, —1-3. 
We h a ve on ly included sources from the IDe Breuck et al] 
|200ll . 120021 ) sample that have meas ured K magn i tudes . 
The combined sample compiled by IWillott et al.l (|2003l ) 
includes 202 narrow- line galaxies from the 3CRR, 6C*, 
6CE and 7CRS catalogues, selected at 151 and 178 MHz. 
Spectroscopic redshifts are available for 193 (96 per cent) of 
those. The 1400 MHz flux densities were obtained from the 
NASA/IPAC Extragalactic Database (NED) by selecting 
the most recently published integrated flux density. In 



4 Flux densities came from ICondon et al. (1998), 
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Figure 9. The observed-frame radio spectral energy distribution 
for NVSS J213238-335318 at z = 2.900. Of the sources in Ta- 
bles \5\ and \E\ this source was fitted to a power-law using the 
largest frequency range, from 151 to 8640 MHz. The resultant 
spectral index and goodness-of-fit statistics are shown. 

several cases only peak flux densities were available and 
25 targets (13 per cent) have no published observations at 
1400 MHz. As the 3CRR/6C/6CE/7CRS combined sample 
has nearly-complete redshift followup, the histograms shown 
in Fig. [10] are not bias ed compa red to the corresponding 
7f-band distributions in lPaper ID . 

It is important to note that while the if-magnitude 
measurements are mostly complete for each of the USS sam- 
ples, the redshifts are incomplete. The percentage of K- 
band detections with spectroscopic redshifts is 32 per cent 
for 6C**, 46 per cent for SU MSS-NVSS, 40 per cent for 
iDe Breuck et all (|200ll . I2002T J and currently 21 per cent for 
the present MRCR-SUMSS sample. 

The choice of targets for follow-up spectroscopy in the 
MRCR-SUMSS sample was essentially determined by the 
telescope being used. On the one hand, the faintest objects 
are the best high- z candidates for followup spectroscopy, but 
the available telescope time will always limit the number of 
faint objects that can be observed. The brightest Tf-band 
objects were observed with the ANU 2.3 m, and those that 
were not detected were then observed on the NTT along 
with fainter targets down to K < 19. After the magnitude 
cut, the targets were then chosen on Right Ascensions avail- 
able. The K > 19 objects were observed at the VLT, where 
we had a slight preference to observe the more compact ob- 
jects, but this was also limited by the accessible Right As- 
censions. To test whether the objects that have spectroscopy 
are representative of the whole sample, we distinguish them 
with different symbols in Fig. [7] which shows largest angu- 
lar size of the sources versus 7f-band magnitude. While the 
distribution of the objects with spectra is similar to that for 
the full sample, among the faintest targets observed on the 
VLT there is an excess of the most compact sources. Larger 
sources could be faint because they are genuinely distant 



lLaine & Peacockl dl98fjf). 


Kellermann 


Paulinv-Toth & Williams! 


119691). IWhite & Beckei 


dl992h. 


Owen & Ledlowl dl997l). 


ICroston et al.l J2005h and 


Becker, White & Hclfandl dl995l). 



or because they are closer but dust obscured. The effect of 
this selection on our redshift distribution would be a small 
bias towards higher-redshift objects. Nevertheless, the me- 
dian 7C-band magnitude for the spectroscopic subset is 18.3 
compared with 18.7 for the full sample (with known QSOs 
removed). This indicates that the spectroscopic sample is 
representative of the full sample, with the slightly brighter 
median simply due to the fact that we have had more NTT 
than VLT observing time. The USS-selected literature sam- 
ples also have median K-band magnitudes for their spectro- 
scopic samples that are ^ 0.2 magnitudes brighter than for 
the full A"-band samples. 

A more stringent test of whether the spectroscopic sam- 
ples are typical of the complete samples is to compare flux 
density distributions. I n Fig. | 10l we compare the flux den- 
sity distributions from IPaper ij (fig. 10) with those from 
the literature. The sources with redshifts from the MRCR- 
SUMSS, SUMSS-NVSS, 6C**, CENSORS and McCarthy 
samples are all consistent with being drawn randomly from 
the flux density distributions of the whole i\-band samples, 
based on Kolmogorov-Smirnov test p values of 0.97, 0.986, 
0.9 10, 0.864 and 0.58 respectivel y. The p value of 0.099 for 
the lDe Breuck et al.1 (|200ll . 120021 ) sample indicates that the 
sources with redshifts from this sample are somewhat less 
representative of the original isT-band population. 

Now that it is clear there are no strong biases in 
the K-band selection or the flux densities of the sources 
that have redshifts, we can investigate the correspond- 
ing redshift distributions. The CENSORS, McCarthy and 
3CRR/6C/6CE/7CRS redshift distributions in Fig.fTOlare 
weighted towards low redshifts, while the USS-selected sam- 
ples have returned higher median redshifts, as expected. The 
different selection frequencies and spectral indices result in 
the variation in flux densities between samples. In IPaper III 
we found that the if-magnitude distribution did not vary as 
a function of flux density. The redshift distributions for the 
USS samples also do not correlate with flux density. On the 
one hand, the SUMSS-NVSS and MRCR-SUMSS samples 
have similar median redshifts but markedly different flux 
densities. On the other hand, the two samples netting the 
highest redshifts have high median flux densities. 

In IPaper il we discussed how the k-correction would 
predict that samples selected at lower radio frequencies 
should result in a higher median redshift. Considering just 
the groups that have been shown to have redshifts repre- 
sentative of the full sample, the MRCR-SUMSS, SUMSS- 
NVSS and 6C** have selection frequencies of 408, 843 and 
151 MHz respectively. For these three samples the median 
redshift does not correlate with selection frequency, sup- 
porting a picture in which the k-correction is not responsi- 
ble for the steep spectral-energy distribution of high-redshift 
sources. 

An alternative explanation for the steep-spectral- 
energy distributio n of radio galaxies was given by 
lAthreva fc Kapahl ljl998h . They suggested that higher- 
density regions at high redshift reduce the hotspot veloc- 
ities through the intergalacti c medium, resulting in steeper 
electron-energy distributions. [Klamer et aD (|2006h then pro- 
posed that the mounting evidence th at HzRGs live in 



over-dense regions or protoclusters (e.g. ICarilli et al.l 


1997|; 


Kurk et al.ll2000bl: IPentericci et al. 2000: Stcidcl et al.1 


2005; 


Kodama et al.ll2007l: Venemans et al.ll2007l) and the fact that 
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their radio spectral energy distributions are straight, implies 
that the z—a correlation is due to the dense environment at 
high redshift. Fig.[TT]shows spectral index versus redshift for 
MRCR-SUMSS and for three USS-selected samples and one 
non-USS-selected sample from the literature. All three USS- 
selected literature samples have a nominal spectral index 
cut-off of a = -1.3, while the MRCR-SUMSS cut-off was 
— 1.0. The non-USS-selected sample is based on the d ata in 
iMcCarthv et al. fl996; 2007. private communication)! , from 
which the sources with redshifts were cross-matched with 
the NVSS catalogue and a spectral index was calculated be- 
tween 408 and 1400 MHz. Fig. 1111 shows an apparent trend 
in spectral index with redshift for the USS-selected sources 
which is not real; it arises because the spectral index cut-off 
in each sample has been breached by objects whose spec- 
tral index changed after sample selection because of revised 
flux densities. Neither the points in the literature samples 
with a < -1.3, nor the points in the MRCR-SUMSS sam- 
ple with a < —1.0 show a significant correlation with red- 
shift. In Fig. 1101 the more stringent spectral index cut-off of 
a < -1.3 of the SUMSS-NVSS sample did not result in a 
significantly higher redshift than did the a < —1.0 selection 
of the MRCR-SUMSS and 6C** samples. However, Fig. [11] 
clearly illustrates that a spectral index cut-off of a = —1 
eliminates the bulk of the low redshift (z < 1) sources. Ap- 
plying a steeper spectral index cut-off will miss many z > 2 
objects. For the objects shown, ~ 54 per cent of the z > 2 
galaxies, or ~ 49 per cent of the z > 2.5 galaxies, would be 
missed by an a — —1.3 cut-off. However, the z-a correla- 
tion is most clearly defined by the lower envelope in the z-a 
plane, which is seen best in the non-USS-selected points. 

If the higher-density environment at high redshift drives 
the steep radio spectra, then it is interesting to note that in 
Fig. 1111 the sources steeper than a < —1.3 are no more 
likely to be at high redshift than at low redshift. However, 
the source numbers are low, so we have also plotted the full 
MRCR-SUMSS K-band sample against spectral index in 
Fig. 1121 While an individual K magnitude is not a direct 
indicator of redshift (see Section 14.1.20 . if the net distri- 
bution of K magnitudes becomes fainter at higher redshift 
then the a versus K plot in Fig. [12] should show a trend 
for sources with a < —1.0. However, there is no evidence 
that the USS sources steepen further with fainter K magni- 
tudes, suggesting there is no advantage in selecting a more 
stringent spectral index cut-off. 

So far, USS-selected samples have been efficient at find- 
ing high-redshift radio galaxies, although non-USS samples 
have also identified many z > 2 obj ects that may well 
have been missed by USS selection (e.g. iBrookes et al.ll2006l . 
2008). While the USS-selection nets higher-redshift samples 
(Fig. I10p . a secondary gain is that USS-selection requires 
substantially less observing to return the same redshift ben- 
efits. For example, CENSORS has 122 sources in a 6 deg 2 
field, while the USS-selected SUMSS-NVSS and MRCR- 
SUMSS samples have 65 and 162 (160 with 4-arcsec-aperture 
magnitudes) 7f-band objects respectively across the sky. 
Considering only the K 19 galaxies which are potentially 
the most distant, the USS-selected samples have > 1.5 times 
the percentage of K 19 objects than the non-USS sam- 
ple. Of the K ^ 19 sources with spectroscopic redshifts, the 
USS-selected samples have a larger fraction of high-redshift 
(z > 2.5) objects, but the numbers are small, especially 
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Figure 11. Spectral index versus spectroscopic red shift for our 
MRCR-SUMSS sample (circles; using a7_ p0 mt from lPaper III ta- 
ble 3 or, for sources without 4 800- and 8640-MHz radio data, we 
used as—point or c?g|^ from |PaperIl| table 2) and thre e sam - 
ples from the literature (stars): |Pe Breuck et al] d200ll . l2002h 
(using «^ 4 3° 5 ° ), SU MSS-NVSS sample (using linear fitted spec- 
tral indices from [K lamer et al] ll2p06t) . or when unavailable, 
Q 843°) and iBornancini et al] J2007I) (alif). The dotted line 
shows the spectral index cut-off of the literature sampl e s (ex- 
cept a small number of sources from iDe Breuck et al] ll200ll . 
l2002t) had a cut-off at a < -1.2), while our MRCR-SUMSS 
cut-off is at a = —1.0. In each sample there are some sources 
which have shallower spectral indices than the cut-offs because 
of subsequent revisions of source flux densities in the cata- 
lo gues. The small dots are from the non-USS-selected sample 
of IMcCarthv et al. (1996; 2007. private com municat ion)! , with a 
spectral index measured between 408 and 1400 MHz (see text for 
details). 



for the SUMSS-NVSS sample. This means that the non- 
USS sample would require at least 1.5 times more K— band 
imaging than the USS samples in order to net a smaller per- 
centage of high redshifts. Therefore, if the aim is principally 
to get a high-redshift sample of objects, then USS selection 
will reduce the amount of imaging required. 



4.5 Spectral line diagnostics 

The bright UV emission lines in HzRGs can be produced 
by one of several ionisation mechanisms. The central AGN 
can photoionise the emission-line region gas, while shock 
ionisation can result from the expanding cocoon associ- 
ated with the radio jets, which shocks the lower tempera- 
ture gas clouds embedded within the higher temperature 
gas surrounding the AGN. In high-redshift galaxies, UV 
emission lines are shifted into the optical band, and can 
be used t o distinguish between sh ock and photoionisation 
processes. IDe Breuck et al] (|2p00b) discuss in detail the dif- 
fer ent models for shock excit ation (based on the models 
of iDopita fe Sutherland! [l996l) and photoionisation (b ased 



on models by Villar-Martm. Tadhunter fe Clark] Il997l and 
I Allen. Dopita fc Tsvetanovll 19981 ) resulting in a set of diag- 
nostic line-ratio plots. In Fig. 1131 we have reproduced those 
models on the plots of the line ratios for our galaxies. There 
are 13 spectra in Fig. [T] which have at least two of the UV 
diagnostic lines CivA1549, HellA1640 and Cm]A1909. For 
the sources which show only two of the three lines, we deter- 
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Figure 10. Distribut i ons o f redshifts and 1 4 00 MHz flux densiti es for the samp les from MR CR-SUMSS, SUMSS-NVS S 
llDe Breuck et al] |2004|. 120061). |Pe Breuck et alj <200lL 120021 ,. 6C** ijCruz et al.l 120061 ). CENSORS ferookes et al] 120061 . 120081) . 
[McCart hT'eTlj"*7lii96T 2007. private communication)! and a combined sample from the 3CRR, 6C*, 6CE and 7CRS catalogues (compiled 
bv lWillott et al.ll2003h IPaper lJ included histo grams of the flux densities of all the sources with i^-band magnitudes, while here we show 
the flux densities of the subset of sour ces with redshifts. The CENSORS redshift histogram includes both the galaxies for which we 
had calculated a 4-arcsec magnitud e in IPaper III as well as t hose that only had 3-arcsec-a perture magnitudes and hence are not in the 
iY-magnitude histogram of IPaper III . Both the McCarthy and lDe Breuck et al.l (2001, 2002) samples only have data shown if the K mag- 
nitude was measured, while sources with a redshift but no K magnitude are not included. This was done because we are comparing the 
redshift distributions of i£"-band-selected galaxies and those which did not have K magnitudes were based on different selection criteria. 
Vertical arrows mark the median values in each plot. QSOs have been removed from the literature samples. The redshift distribution for 
the MRCR-SUMSS sample shows both the radio galaxies and QSOs (dotted line) and just the radio galaxies (solid line), but the median 
calculation did not include the QSOs, in order to be directly comparable, to the literature samples. 
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K— band magnitude (4" aperture) 

Figure 12. Spectral index versus _R"-band magnitude in a 4-arcsec 
aperture for our MRCR-SUMS S sam ple. The spectral indices are 
the 7-point fitted values from lPaper III table 3 for sources that have 
4800- and 8640-MHz radio data, otherwise we used the 5-point 
fit ted value s or, when unavailable, the Ogljj values as defined 
in iPaper III (table 2). Our spectral index cut-off was defined us- 
ing 408—843 MHz, and therefore some of our sources turned out 
to be flatter when fitted over a large frequency range, and some 
had the 843-MHz flux densities revis ed by th e new SUMSS cat- 
alogue (version 2.0), as discussed in Pa per j The sources with 
_R"-band magnitude limits are marked by a circle rather than a 
right-pointing arrow, for clarity. For sources steeper than —1.0, 
there is no correlation between spectral index and K magnitude. 



mined 3a upp e r limi ts to the flux from the missing line. Fol- 
lowing iHobbsl | 1984t ). the limiting flux was calculated from 
the product of three times the rms fluctuation in the con- 
tinuum and a line width Av. For He II we set Aw equal to 
the width of the C IV line (since they are both permitted 
lines with similar ionisation potentials), while for C ill] we 
set Av equal to the instrumental resolution, since this line is 
typically narrow in radio galaxies. All of the galaxies plotted 
have z > 1.8. 

The diagnostic diagrams show that all of the galax- 
ies are consistent with photoionisation as the primary ex- 
citation process. There are no galaxies which can be ex- 
plained by shocks alone or shock plus precursor excita- 
tion, but there is one with a limiting arrow which may 
extend to the shock models. The majority of our galaxies 
lie between the lines representing photoionisation models 
where the incident photoionising continuum has a power 
law spectral index of a > —1.5 and ionisation parameters 
of —2.5 < log 10 (f7) < —1.5. These result s agree with the 
much la r ger sample of HzR Gs presented in lDe Breuck et aU 
l|2000bl ). iMcCarthvl (|l993l ) produced a composite spectrum 
for 0.1 < z < 3 radio galaxies, which also demonstrated the 
dominance of photoionisation in radio galaxies. 

The QSOs in Fig. \T3\ typically lie above the galaxies 
in each plot due to a comparatively smaller Hell line flux. 
The co mposite QSO spectrum fr om the Sloan Digital Sky 
Survey (|Vanden Berk et al.ll2CX)iT ) has log(C m]/Heii) = 1.5 
and log(C iv/He n) = 1.7, which are very much higher than 
the QSOs in our sample. The QSOs in our sample have ratios 
much closer to that of the radio galaxies. Compact Steep- 
Spectrum (CSS) radio QSOs from lBaker fc Hunsteadl (|l995l ) 
similarly showed stronger Hell than typical QSOs. 



In [Paper III we found that some sources in our sam- 
ple show alignment between the radio structures and either 
the major axis of the fitted Jf-band ellipsoid, or the linear 
extended Tf-band emission. It was proposed that this align- 
ment could be due to jet-induced star formation or alterna- 
tively, scattering of the nuclear emission along paths carved 
out by the jets through the gas. NVSS J210814-350823, 
NVSS J213637-340318 and NVSS J223101-353227 have an 
ellipsoid major ax is aligned within 10 degrees of the radio 
axis (see IPaper IS ) , and are included on the emission line 
diagnostic plots in Fig. 1131 All three have line ratios that 
do not match the shock models. On the one hand, we may 
expect shock excitation in regions of jet-induced star forma- 
tion, but on the other hand, our spectra were extracted in 
~ 1-arcsec apertures, which will be dominated by the nu- 
clear source. Therefore we can not rule out jet-induced star 
formation on the basis of the line ratios alone. 

The spectra of HzRGs are typically dominated by the 
Lya line in the UV spectrum. Lya halos can ext end to hun- 
dreds of kiloparsecs around the radio source (|Kurk et al.l 
l2000al ; IVenemans et al J 120021 ; iReuland et~al1l2007h . Hi gas 
surrounding the active nucleus, or along the line of sight, 
can abs orb Lya, cutting off the blue wing of the spec- 
tral line. Ivan Oiik et al.l l|l997h and lDe Breuck et all (|2000bl. 
fig. 11) have shown that high-redshift galaxies with absorp- 
tion in the Lya line tend to have compact radio structures 
(< 50kpc). While the resolution of our spectra is insuffi- 
cient to examine the Lya profile in detail, there are two 
quasars in the MRCR-SUMSS sample with obvious ab- 
sorption in the blue wing of the Lya line. The absorption 
profiles are shown in Fig. M NVSS J105917-303658 and 
NVSS J144206— 345115 have compact linear sizes of 7.6 and 
7.2 kpc, and have broad lines. The small apparent size may 
be due to alignment of the jets to the line-of-sight, or the jets 
may be intrinsically compact. In the latter case, these two 
source s are consistent with t he find ings of Ivan Oiik et ail 
< |l997f ) and |Pe Breuck et all (|2000tJl that stronger Hi ab- 
sorption of the Lya line may be due to a denser medium 
surrounding and confining the radio source. An alternative 
model involves a low-density medium with an absorbing 
halo of low-metall icity gas which has not yet been impacted 
by th e radio jets l|Binette. Kurk. Villar-Martm. Rottgering 
l200Cf ). Evidence for denser environments around these and 
other sources will be discussed further in a subsequent paper 
based on Faraday rotation measures and galaxy overdensi- 
ties in the surrounding fields. 

Dust in elliptical galaxies is typically associated with 
recent mergers. If HzRGs are forming through hierarchi- 
cal mergers, which may be the trigger for the nuclear 
engine and hence the radio jets, then we might expect 
the most compact, newly-triggered sources to show dust 
extinction. We now consider the evidence for dust from 
spectral lines. While Lya is generally the strongest UV 
emission line in HzRG spectra, it is not a good mea- 
sure of dust as it is more likely to be be weak relative 
to the other UV lines as a result of viewing angle ef- 
fects rather than dust (|Villar-Martm. Binette fc Fosburvl 
1996) . Based on the dust extinct ion vectors calculated by 
Allen. Dopita fc Tsvetanovl (Il998l . fig. 4), the C m]A1909 line 
is more affected by dust than C IVA1549. Therefore, in spec- 
tra such as that for NVSS J215226-341606, the marginal 
detection of Cm]A1909 may be indicative of dust obscura- 
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log (CIII] A,1909/HeH M640) log (CIV A.1549/CIII] M909) 




log (CIV A.1549/CIII] A.1909) 

Figure 13. UV spectral line diagnostic plots for the lines CivA1549, Cin]A1909 and HenA1640. The galaxies in our sample that have 
at least two of these lines have been plotted (dots). In some cases an upper limit has been used for one of the lines. Radio galaxies are 
marked by dots, while th e QSOs have square symbols. We have also plotted photoionisation and shock models reproduced directly from 
|Pe Breuck et al. (2000b); the details of the model parameters are discussed in that paper. The light grey shaded grid represents shock 
models while the unshaded grid (which partially overlaps the shock models) shows the shock+precursor models. The shock velocities 
start at 150 kms" 1 at the top right of the curves, and increase along the curves to 500 kms" 1 at the end. The four curves in each 
grid represent four values of the magnetic parameter < B/y/n < 4/iGcm ' 2 . Shock models for 500-1000 km s" 1 with n = 1.0 and 
B = 3.23 from lAllen et al.l (2008) are marked by a dashed line. The dark grey shaded grids represent photoionisation models with four 
photo-ionisation sequences, with values of the ionisation parameter log 10 (£/) every 0.5 dex. The top and bottom shaded photoionisation 
grids have two models each, with power law spectral indices of a = — 1 and a = —1.5 respectively. The left or bottom edge of the 
photoionisation grids are models with hydr ogen density n = 100 cm" 3 and the top or right edges have n = 1000 cm" 3 . The single line 
with crosses along it represents the model of lBinette. Wilson &: Storchi-Bcrgmann ( 1996). This model assumes that the photoionised light 
has gone through both matter-bound and ionisation-bound clouds in which the ratio of the solid angle subtended by the matter-bounded 
to ionisation-bounded clouds is varied from 0.01 to 100 alo ng the line with tick marks e v ery .02 dex. An arrow represents 10 magnitudes 
of extinction based on the Milky Way extinction law from lAllen. Dopita &: Tsvetanovj [|l998h . Two QSOs (NVSS J152123- 375708 and 
NVSS J152435-352623) have upper limits on the flux of their HellA1640 lines which give Civ/Hell> 1 and for NVSS J152435-352623 
Cni]/Hcn> 1 and therefore these values are outside the region plotted. 
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Figure 14. Line profiles of Lya and CivA1549 on a common 
velocity scale for the two sources which show self-absorption, 
NVSS J105917-303658 and NVSS J144206-345115. The posi- 
tion of the common absorption is marked with the dashed line. 
The velocity origin was determined by the redshift in Table [2] 



tion, although there are insufficient other lines to confirm 
this. The only other galaxy in which Clll]A1909 was not 
detected when HeiiA1640 and/or CivA1549 were detected 
is NVSS J144932-385657. In this case, based on an up- 
per l imit on the Ciii1A1909 line , the d iagnostic diagrams 
from lAllen. Dopita fc Tsvetanovl (|l998l . fig. 4) can explain 
the lack of Cm]A1909 with either photoionisation or shock 
models, without the need for large dust extinction. 



5 CONCLUSIONS 

Spectroscopic results and analysis of the radio, is'-band and 
redshift data have been presented for a large Southern hemi- 
sphere USS-selected HzRG survey, the MRCR-SUMSS sam- 
ple. Our main findings are: 

(i) Based on 175 images and 52 spectra, 164 radio 
sources have A'-band or SuperCOSMOS identifications with 
36 confirmed galaxy redshifts. The highest new redshift 
found is 3.26, but our spectroscopic followup is so far only 
30 per cent complete. A comparison of the USS samples 
with the non-USS-selected CENSORS and McCarthy sam- 
ples shows that the fraction of z < 1 galaxies is signifi- 
cantly reduced by USS-selection, as expected. However, for 
the USS-selected samples there is no correlation between the 
median redshift and either median flux density or selection 
frequency. 

(ii) Detailed analysis of the K-z distribution includes 
a fit to 4-arcsec-aperture K magnitudes for several sam- 
ples from the literature. Th i s fit is co mpared with the 
64-kp c fits of IWillott et all (|2003T ) and |Pe Breuck et all 
(l2004h. The MRCR-SUMSS sample agrees with the 
iRocca-Vo lmcrangc ct al. | l|2004h models for > lO n M galax- 
ies for all but two sources. We find no evidence for an in- 
crease in dispersion in the K-z plot at z > 2 and conclude 
that these galaxies have undergone passive stellar evolution 
since z > 3. 

(iii) To test for curvature in the radio-spectral-energy 



distribution at frequencies lower than our selection fre- 
quency, both the MRCR-SUMSS and the SUMSS-NVSS 
samples were cross-matched with the 74-MHz VLSS cata- 
logue and 151-MHz MRT catalogues. 75 per cent (15/20) 
of the sources with if-band magnitudes have radio-spectral- 
energy distributions that can be fitted by a single power-law. 

(iv) The k-correction has frequently been assumed to 
account for the z — a correlation and hence for the success 
of USS-selection in finding high-redshift-radio galaxies. The 
predominance of straight spectrum objects over frequency 
ranges up to 74-8640 MHz and the lack of correlation with 
selection frequency are both evidence that an entirely differ- 
ent mechanism is responsible for the steep spectra of high 
redshift radio galaxies. Suggestions by iKlamer et al l ( 2006) 
that the steep spectrum is due to the higher-density environ- 
ments at high redshift will be tested for the MRCR-SUMSS 
sample in a subsequent paper. 

(v) The spectral line ratios indicate that photoionisa- 
tion is more dominant than shocks as the ionisation mecha- 
nism at the centre of these galaxies. QSOs in our sample have 
line ratios that are much closer to that of the radio galaxies 
than typical QSO ratios due to a stronger HenA1640 line, 
which is also seen in CSS sources. We find no evidence for 
significant dust extinction around the nucleus of any of our 
sources based on the line ratios. 

Our method has been successful in identifying twelve 
2 < 2 < 3.5 galaxies which are at the epoch when the 
star formation rate density of the Universe was at a peak. 
They form an important sample for investigating the envi- 
ronments and evolution of massive galaxies. By the time the 
spectroscopic followup is complete, we expect more than 40 
galaxies in this redshift range. In an upcoming paper, we 
will present clustering and evolution studies based on rota- 
tion measures and source overdensities in our images. 

This paper has been typeset from a TffjX/ ]M£X file prepared 
by the author. 
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